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Coherent  Thomson  scattering  of  a  picosecond  probe  laser  was  used  to  measure  the  time  evolution  of  plasma  wakefields  produced  by  a 
high  intensity  laser  pulse  (4  x  1018  W/cm2)  in  an  underdense  plasma  (ne  m  1019  cm'3)  in  the  Self-Modulated  Laser  Wakefield 
Accelerator  configuration.  Plasma  wakefields  are  observed  to  last  for  approximately  5  psec.  Ion  acoustic  waves,  generated  by  the 
decay  of  the  plasma  Wakefield,  were  measured  to  have  a  duration  of  approximately  50  psec.  These  observations  are  consistent  with 
the  decay  of  large  amplitude  plasma  waves  due  to  the  modulational  instability  and  the  enhancement  of  the  scattered  signal  from  plasma 
channel  formation. 
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TEMPORAL  EVOLUTION  OF  SELF-MODULATED  LASER  WAKEFIELDS 
MEASURED  BY  COHERENT  THOMSON  SCATTERING 


Significant  progress  has  been  made  in  recent  years  |1|  in  the  use  oi  laser-produced 
plasma  as  a  medium  for  accelerating  electrons  to  high  energies.  In  the  Laser  Wakefield 
Accelerator  (LWFA)  |2|,  a  high  intensity  laser  pulse  is  focused  into  an  underdense  plasma  with  a 
pulse  duration,  Tl  ,  similar  to  the  electron  plasma  period  (i.e.,  Tl  ~  27t/tOpe  where  C0pC  is  the 
electron  plasma  frequency).  Large  amplitude  plasma  waves  (wakefields)  are  generated  with 
strong  longitudinal  electric  fields  and  relativistic  phase  velocities  which  are  capable  of 
accelerating  injected  electrons.  For  laser  powers,  P  ,  approaching  or  exceeding  the  relativistic 
self-focusing  threshold  (i.e.,  Pc  =  17  (copc/030^  GW;  where  co0  is  the  laser  frequency),  it  is  not 
necessary  to  match  the  pulse  duration  to  the  plasma  period.  At  laser  powers  for  which  P  >  Pc  and 
tL  >  27t/C0pe,  the  laser  envelope  undergoes  an  instability  and  becomes  "self-modulated"  at  the 
plasma  frequency  [3].  This  effect  resonantly  enhances  the  creation  of  wakefields  and  allows  use 
of  higher  electron  densities  and  generates  stronger  accelerating  fields  [4],  Recent  experiments  in 
this  Self-Modulated  Laser  Wakefield  Accelerator  (SM-LWFA)  regime  have  measured  the 
production  of  high  energy  electrons  where  the  source  of  accelerated  particles  was  either 
background  electrons  from  the  target  plasma  [5]  or  electrons  injected  into  the  interaction  region 
from  an  adjacent  laser-produced  plasma  [6].  Direct  observations  of  wakefields  in  the 
conventional  LWFA  configuration  were  recently  reported  using  interferometric  techniques,  in 
which  the  spatial  and  temporal  waveforms  of  the  wakefield  behind  the  pump  laser  pulse  were 
measured  [7]. 

This  letter  will  discuss  the  first  experiments  to  use  coherent  Thomson  scattering  (CTS)  of 
a  picosecond  probe  laser  pulse  to  measure  the  temporal  behavior  of  self-modulated  wakefields  in 
the  SM-LWFA.  We  confirmed  the  generation  of  wakefields  with  phase  velocities  close  to  the 
speed  of  light.  We  observed  the  turbulent  decay  of  the  wakefields  and  the  transfer  of  energy  to 
ion  acoustic  waves,  possibly  through  a  parametric  process  such  as  the  modulational  instability. 
The  measured  decay  time  of  these  waves  is  consistent  with  growth  rates  of  the  modulational 
instability  for  our  experimental  conditions.  The  formation  of  a  plasma  density  channel  behind 
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the  pump  laser  pulse,  which  can  guide  the  probe  laser  pulse,  could  also  enhance  the  Thomson 
scattered  signal. 

The  Tabletop  Terawatt  laser  at  the  Naval  Research  Laboratory  is  a  Ti:Sapphire/Nd:Glass 
system  (X  =  1.054  pm)  which  uses  chirped  pulse  amplification  (8]  to  achieve  the  high  powers 
necessary  for  experiments  in  the  self-modulated  regime.  In  these  experiments,  the  pump  laser 
pulse  had  a  typical  duration  of  400  fsec  and  an  energy  of  800  mJ  (P  =  2  TW).  The  contrast  (i.e., 
the  ratio  of  the  prepulse  intensity  to  that  at  the  peak  of  the  laser  pulse)  was  measured  to  be  10‘6. 
When  the  beam  was  focused  at  1  TW  by  an  off-axis  parabolic  mirror  (f/3),  the  measured  focal 
spot  diameter  in  vacuum  was  8.5  pm.  This  implies  a  peak  focused  intensity  in  our  experiments 
of  7  x  1018  W/cm2.  For  a  plasma  density  of  n0  ~  1()19  cnv3,  the  critical  power  is  Pc  ~  1.8  TW, 
and  the  plasma  wavelength  is  Xpe  =  2TCc/(Ope  ~  10  pm.  Hence,  P  >  Pc  and  ctp  ~  12  Xpe>  which  are 
necessary  for  operation  in  the  SM-LWFA  regime. 

In  CTS,  waves  in  a  plasma  will  enhance  Thomson  scattered  light  at  particular 
wavelengths  and  directions  due  to  interference  effects.  Coherently  scattered  light  has  a  wave 
vector  ksc  and  a  frequency  cosc  that  satisfy  the  Bragg  scattering  conditions  of  frequency  and 
wavenumber  matching.  For  electron  plasma  waves  with  G)pe  and  kpe,  the  conditions  require  (%c 
=  0)i  ±  0)pe  and  ksc  =  k  j  ±  kpe,  where  coi  is  the  frequency  and  k]  is  the  wave  vector  of  the  probe 
laser.  Plasma  waves  capable  of  accelerating  electrons  to  high  energy  (as  produced  by  a  LWFA) 
propagate  collinearly  with  the  laser  pulse  and  have  relativistic  phase  velocities,  v,),  ~  c.  CTS  has 
also  been  used  to  measure  similar  relativistic  (v^  ~  c)  plasma  waves  in  the  laser  beat-wave 
accelerator  [9,  10].  For  correct  matching  of  k  vectors,  both  the  probe  and  the  Thomson  scattered 
light  must  propagate  in  the  same  direction  as  the  relativistic  plasma  wave  (i.e.,  along  the  pump 
axis  of  propagation).  For  coherent  scattering  to  be  observable,  the  condition  (k  Xd)  <  1  must  be 
satisfied,  where  X&  is  the  plasma  Debye  length  and  k  is  the  wave  number  of  the  plasma  waves  to 
be  measured  [1 1]  -  otherwise  the  waves  undergo  severe  Landau  damping.  In  these  experiments, 
we  estimate  ~  0.01  pm  and  k  ~  o)pe/c  ~  0.6  pnr1. 
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The  setup  for  our  experiments  is  shown  in  Figure  1.  Approximately  10  %  of  the  main 
beam  was  split  off  by  a  pellicle  for  use  as  a  probe  pulse  and  was  subsequently  frequency  doubled 
(to  527  nm)  by  a  1  cm  thick  KD’P  crystal.  For  800  mJ  pump  pulses  the  probe  pulse  energy  was 
typically  10  mJ.  The  maximum  broadening  of  the  probe  pulselength  from  passage  through  the 
doubling  crystal  will  result  in  a  pulse  duration,  Tl,p  <  1  psec.  This  provides  the  intrinsic 
temporal  resolution  for  our  experiment.  Since  ctl.p  »  ^-pe  in  these  experiments,  the  probe 
pulse  efficiently  sampled  the  plasma  waves  it  measured.  This  is  in  contrast  to  optical 
interferometric  techniques  [7]  which  require  much  shorter  probe  pulses  such  that  ctL.p  <  ^-pe- 

The  probe  pulse  traverses  a  delay  line  and  then  is  recombined  with  the  main  beam  before 
being  focused  into  the  plasma  created  by  the  pump  pulse.  This  delay  line  can  be  adjusted  to  vary 
the  timing  between  the  two  pulses.  The  paths  of  the  two  beams  were  overlapped  by  alignment  to 
a  series  of  apertures  prior  to  the  focusing  optics  (f/3  off-axis  parabolic  mirror)  as  well  as  by 
observation  of  scattered  light  in  the  focal  region.  The  target  was  a  gas  jet  of  either  helium  or 
hydrogen  and  we  measured  the  electron  density  by  recording  the  wavelength  shift  of  stimulated 
Raman  backscattered  light  from  interactions  at  lower  intensities  (lO1^  W/cm^).  The  electron 
density  used  in  these  experiments  was  in  the  range  of  0.5  -  1.2  x  lO1^  cm*3  with  ~  25  %  of 
accuracy.  The  laser  breakdown  is  monitored  with  a  CCD  camera  at  90°  to  the  propagation  axis 
of  the  pump  laser. 

Coarse  temporal  synchronization  of  pump  and  probe  pulses  was  achieved  using  a  streak 
camera  with  better  than  4  psec  resolution.  Fine  temporal  synchronization  was  attained  by 
measuring  the  probe  frequency  after  passage  through  the  interaction  region.  If  the  probe  pulse 
travels  ahead  of  the  pump,  it  ionizes  the  gas  and  generates  a  blue  shift  in  the  spectrum  [12]. 
Otherwise,  the  frequency  is  unshifted.  This  effect  is  easily  observed  and  provides  a  measurement 
of  the  temporal  overlap  to  within  1  psec. 

These  experiments  used  two  different  CTS  configurations.  The  first  configuration  had 
the  pump,  probe,  and  scattered  light  propagating  collinearly  and  unscattered  light  from  the  pump 
beam  was  attenuated  by  an  infrared  absorption  filter  after  the  interaction  region.  The  majority  of 
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the  527  nm  probe  light  was  not  scattered  and  was  blocked  before  the  slit  of  the  spectrometer 
(0.25  m  Czerny-Turner)  by  a  notch  filter  in  the  beam  path.  However,  the  Thomson  scattered 
electron  plasma  satellites  of  the  probe  light,  which  are  shifted  by  0)pC  (AX-siokcs  ~  23  nm,  AXgnii  - 
Stokes  ~  21  nm  for  a  plasma  density  of  7  x  1018  cm-3  as  in  fig.  2),  are  positioned  beyond  the 
edges  of  the  absorption  band  of  the  notch  filter  and  hence  will  be  detected  by  the  spectrometer. 
The  principal  result  in  this  configuration  was  the  observation  of  these  plasma  satellites  for  about 
5  psec  after  passage  of  the  pump  laser  (see  Figure  2,  which  shows  a  FWHM  of  ~2  psec).  As 
shown  in  the  insert  in  Fig.  2,  the  anti-Stokes  line  was  typically  more  intense  than  the  Stokes  line, 
indicating  that  electron  plasma  waves  in  the  wakefield  preferentially  propagate  in  the  forward 
direction  [13].  These  measurements  confirmed  the  generation  of  wakefields  with  v<()  ~  c  [14]. 

The  ratio  of  the  total  scattered  signal  intensity  to  the  probe  intensity  is  given  by  [15], 
Is  /  1 1  ~  2.5  (8n  Lz  /  ncr  A.])2,  where  ncr  is  the  critical  density  for  the  probe  wavelength  (4  x  1021 
cm'3).  A.]  is  the  probe  wavelength,  Lz  is  probe-plasma  interaction  length  (the  confocal  parameter 
of  the  probe  laser,  ~  200  pm)  and  5n  is  the  density  perturbation  of  the  wave.  We  estimate  the 
ratio  (Ig/Ii)  to  be  10'3±1  for  scattering  into  the  electron  plasma  satellites.  This  implies  that  the 
amplitude  of  plasma  waves  encountered  by  the  probe  pulse  (8n/n0i  where  n0  is  the  ambient 
density)  was  up  to  10' This  is  probably  an  underestimate  of  the  actual  amplitude  of  the 
wakefields  produced  in  these  experiments  since  slight  misalignments  in  our  setup  can  result  in  a 
substantial  reduction  of  signal. 

In  the  second  CTS  configuration,  the  scattered  light  was  observed  at  a  convenient  forward 
scattering  angle  of  30°  ±  5°  from  the  pump-probe  axis  by  relocating  the  collection  optics  shown 
in  fig.  1.  The  wakefields  generated  in  these  experiments  are  highly  three  dimensional  due  to  the 
tight  focusing  of  the  pump  laser.  This  causes  a  strong  radial  profile  of  the  wakefield  which 
implies  that  the  scattering  plasma  waves  have  significant  kj_  components  (~  27t/ro  ,  ro  being  the 
pump  laser  focal  spot  radius)  and  lower  phase  velocities.  At  30°,  the  phase  velocity  can  vary 
from  0.06  c  to  0.49  c  from  consideration  of  the  f-numbers  of  the  alignment  and  detection  optics. 
These  lower  phase  velocity  plasma  waves  couple  more  readily  with  the  probe  laser  to  scatter  in 
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the  near  forward  directions.  In  fact,  a  10°  tilt  between  the  pump  and  probe  lasers  by  displacing 
the  probe  laser  on  the  focusing  parabolic  mirror  (while  maintaining  the  30°  angle  to  the  pump  for 
the  scattered  light)  gave  a  better  signal  to  noise  ratio  without  changes  to  the  general  behavior  of 
the  scattered  signal. 

In  Fig.  3,  the  Stokes  satellites  at  delay  =  -10  ps  were  generated  by  the  pump  laser  from 
its  interaction  with  plasma  density  gradients  due  to  filamentation  and  charge  displacement  effects 
[16].  When  the  probe  laser  was  injected  at  delay  =  +1  ps,  there  is  a  marked  increase  in  the 
emission  of  the  first  order  electron  satellites.  Electron  satellite  emission  lasted  for  about  5-7 
psec  in  both  hydrogen  and  helium,  similar  to  that  observed  in  the  directly  forward  direction,  as 
shown  by  the  closed  circles  in  Fig.  4. 

Fig.  3  also  shows  a  dramatic  rise  in  the  central  peak  of  the  Thomson  scattered  spectrum 
which  correlates  to  the  decline  of  the  electron  plasma  sidebands.  We  believe  this  central  peak  is 
due  to  coherent  Thomson  scattering  from  ion  acoustic  waves  since  both  a  strong  dependence  on 
the  polarization  of  the  probe  beam  and  a  dependence  of  its  intensity  on  ne2  were  observed.  The 
open  circles  in  Fig.  4  show  that  the  amplitude  of  the  central  scattered  peak  reached  a  maximum  at 
At  =  25  psec  after  the  main  pulse  for  hydrogen  and  at  At  =  40  psec  for  helium.  We  believe  this  is 
due  to  the  turbulent  decay  of  the  forward  going  large  amplitude  plasma  waves  into  ion  acoustic 
waves  through  mechanisms  such  as  the  modulational  instability. 

Significant  theoretical  work  has  been  performed  on  the  decay  of  large  amplitude 
Langmuir  waves  via  parametric  coupling  to  ion  acoustic  waves  [17,  1 8 j .  In  particular,  for 
moderate  temperatures  (Tc<  1  keV),  a  Langmuir  wave  is  subject  to  the  modulation  instability 
(MI)  [17],  in  which  a  plasma  wave  (cop,  kp)  decays  into  a  low-frequency  ion  wave  (G)j,  kj)  and 
two  daughter  electron  plasma  waves  (cop±  (Oj,  kp  ±  kj),  where  typically  kj  »  kp.  Ion  waves 
generated  by  MI  have  been  observed  coincident  with  plasma  waves  in  experiments  on  the  plasma 
beat  wave  accelerator  [11].  In  these  experiments,  the  plasma  wave  is  simultaneously  driven  by 
the  long  duration  laser  pulses  (~  150  ps)  while  decaying  via  MI.  In  SM-LWFA  experiments,  on 
the  other  hand,  the  plasma  waves  are  first  generated  very  rapidly  (<1  ps)  and  then  decay  into  ion 
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waves  in  the  region  behind  the  laser  pulse.  Collisional  damping  of  the  wakefield  can  be 
neglected,  since  for  electrons  oscillating  in  a  8n/n0  ~  0.1  plasma  wave,  the  electron-ion  collision 

time  is  long  (typically  ~  300  psec). 

The  linear  growth  rate  [18]  of  Ml  depends  on  the  parameter  vl/vic  =  23(5n/n0)Tc 
where  Te  is  the  electron  temperature  in  keV,  vL  is  the  electron  quiver  velocity  in  the  Langmuir 
wave  and  v^  is  the  electron  thermal  velocity.  Ml  occurs  in  the  strong  field  limit  when  v^Aie  > 
(C0pe/t0pi)1/3,  where  cope  and  copi  are  the  plasma  frequencies  for  the  electrons  and  ions 
respectively.  Initially,  if  the  plasma  wave  amplitude  is  large  (8n/n0  ~  0.1)  and  for  Te  ~  10  eV, 
the  strong  field  condition  is  satisfied  such  that  vl/viC  ~  23  and  ((Ope/cOpj)*/3  ~  3.  The  maximum 
strong  field  growth  rate  [18]  is  approximately  given  by  the  ion  plasma  period  (  Ysf  ~  (Dpi )  and  the 
corresponding  e-folding  time  xsf  =  1/Ysf  is  short,  about  0.3  psec  for  H2  (2jt/(Dpj  =1.8  psec)  and 
about  0.4  psec  for  He  (2rc/(0pj  =  2.5  psec). 

As  the  initial  wakefield  decays,  however,  8n/n0  decreases  and  MI  transitions  [18]  to  the 
slower  growth  of  the  weak  field  limit,  vl/vic<  ((Dpe/copj)1/3  ,  with  a  maximum  weak  field  growth 
rate  of  Ywf  ~  (Dpj(vp/v^).  Notice  that  as  8n/no  (~  vj^/v ^)  decreases,  so  does  the  growth  rate. 
In  addition,  the  decay  of  the  wakefield  behind  the  laser  pulse  is  likely  to  be  accompanied  by 
electron  plasma  heating,  which  will  also  decrease  the  growth  of  MI.  The  daughter  plasma  waves 
may  also  decay  and  hence  can  lead  to  a  cascade  of  the  wakefield  energy  into  smaller  amplitude 
daughter  waves  and  drive  ion  waves  over  longer  time  scales.  For  example,  Sn/n^  ~  10  and  Tg 
=100  eV  give  xwf  =  1/Ywf  ~  8  psec  for  H2  (1 1  psec  for  He)  which  agree  reasonably  well  with  the 
time  scales  shown  in  Fig.  4  (a)  and  (b).  Another  possible  source  for  ion  waves  is  stimulated 
Brillouin  backscattering  (SBS)  within  the  laser  pulse  [19].  However,  since  the  laser  pulse  is 
short,  the  number  of  e-foldings  for  SBS  will  be  small  ( ~  5  for  typical  experimental  parameters). 

The  temporal  behavior  of  the  central  scattered  peak  shown  in  Fig.  4  may  also  be  affected 
by  the  formation  of  a  plasma  channel.  Detailed  measurements  of  plasma  channel  formation  in 
these  experiments  are  discussed  elsewhere  [20].  As  a  brief  account,  Thomson  scattered 
emission  from  the  pump  laser  was  imaged  onto  the  90°  CCD  camera  shown  in  fig.  1.  An 
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extended  propagation  channel  of  about  2.5  mm  (about  the  width  of  the  gas  jet)  was  observed  as 
the  pump  laser  power  was  increased  to  ~2  TW.  When  the  probe  was  injected  along  with  the 
pump  (Pprobc  «  Pc).  Thomson  scattered  emission  from  the  probe  pulse  at  527  nm  was  observed 
throughout  the  region  where  channeling  of  the  high  intensity  pump  pulse  occurred.  Its  intensity 
increased  to  a  maximum  for  probe  delays  of  ~  15  psec  and  subsequently  decayed  over  about  50 
psec  [20],  This  indicates  that  a  long  lived  waveguide  structure  [21]  may  have  been  created 
which  guided  the  probe  pulse. 

The  plasma  waveguide  is  probably  formed  by  the  following  mechanism.  As  the  pump 
pulse  propagates,  its  ponderomotive  force  radially  expels  electrons.  This  creates  a  strong  space 
charge  field  which  initiates  an  outward  radial  motion  of  the  ions.  After  the  passage  of  the  pump 
pulse,  the  ions  continue  to  drift  radially.  Estimates  indicate  that  the  maximum  depth  of  the 
plasma  channel  An  is  large,  An/n0  >  0.5,  as  is  necessary  to  guide  a  tightly  focused  (ro  ~  5  |im) 
probe  laser  pulse.  The  time  scale  for  the  formation  of  a  plasma  channel  can  be  shown  [20]  to  be 
on  the  order  of  tCh  ~  2'3^2  ro/Cs  ,  where  Cs  =  (ZTe/nij)1''2  is  the  ion  sound  velocity,  Z  is  the  ionic 
charge,  and  m;  is  the  ion  mass.  For  example,  for  He  with  Te  =  100  eV  and  rq  =  5  |im,  it  is  found 
that  tCh  ~  26  psec  (18  psec  for  hydrogen).  Since  the  plasma  channel  guides  the  probe  laser  pulse 
for  longer  propagation  distances  Lz  and  Is  ~  Lz2,  the  CTS  signal  will  be  enhanced  as  the 
channel  is  being  formed.  There  is  reasonable  agreement  between  the  time  scales  for  plasma 
channel  formation  and  the  evolution  of  the  CTS  central  peak  signal  shown  in  Fig.  4. 

In  conclusion,  these  experiments  have  resulted  in  the  first  direct  measurement  of  plasma 
waves  generated  by  a  LWFA  in  the  self-modulated  regime  (n0  ~  lO^cnv3)  using  coherent 
Thomson  scattering  of  a  short  (<  1  psec)  frequency-doubled  probe  laser  pulse.  By  delaying  the 
probe  pulse  with  respect  to  the  pump  pulse,  the  time  evolution  of  the  plasma  wave  was 
measured.  The  lifetime  of  the  wakefield  was  observed  to  be  about  5  picoseconds.  The 
generation  of  ion  waves  was  also  measured  over  longer  timescales,  and  is  consistent  with  the 
theory  of  parametric  decay  of  electron  plasma  waves  via  the  modulation  instability  and  the 
enhancement  of  the  signal  due  to  possible  plasma  channel  formation. 
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[Fig.  1] 


Schematic  of  Thomson  scattering  experimental  setup. 
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[Fig.  2]  Evolution  of  electron  satellites  from  zero-degree  Thomson  scattering.  The 
background  level  is  shown  by  the  dashed  line.  Data  are  not  shown  for  pump- 
probe  delay  times  less  than  zero  since  blueshifting  of  the  probe  laser  spectrum 
saturates  the  detector  for  early  times.  Insert  is  a  typical  spectrum  (Stokes  at  506 
nm,  anti-Stokes  at  550  nm). 
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[Fig.  3]  Evolution  of  Thomson  scattered  spectrum  at  30°  scattering  angle  w.r.t.  pump  laser 
axis  (Pump  and  probe  laser  are  collinear)  for  different  delay  settings  between 
pump  and  probe  pulses.  Negative  delay  implies  emission  from  pump  laser  alone. 
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